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Abstract 

The two major disturbances in the heliosphere during the present sunspot 

cycle, the event of June - August, 1982 and the event of April - June, 1978, 

are simulated by the method developed by Hakamada and Akasofu (1982). 

Specifically, we attempt to simulate effects of six major flares from three 

active regions in June and July, 1982 and April and May, 1978. A comparison 

of the results with the solar wind observations at Pioneer 12 (~ 0.8 au) , - 

• 

ISEE-3 (~ 1 au) , Pioneer 11 (~ 7-13 au) and Pioneer 10 (~ 16-28 au) suggests 
that some major flares occurred behind the disk of the sun during the two 
periods* Our method provides qualitatively some information as to how such a 
series of intense solar flares can greatly disturb both the inner and outer 
heliospheres. A long lasting effect on cosmic rays is discussed in 


conjunction with the disturbed heliosphere. 
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1. Introduc tlon • 

, The months of June and July, 1982 and of April and May, 1978 were two of 
the* most active periods of the sun during the present solar cycle. There were 
at least three active regions during the first period and one region which 
produced a large number of flares, causing major disturbances in the 
' heliosphere. For the first event, several shock waves and/or cosmic ray 
decreases were observed by deep space probes, Pioneer 12 (Venus orbiter) at 
about 0.7 au, Pioneer 11 at about 12 au and Pioneer 10 at about 28 au, as well 
as by the ISEE-3 at the libration point. A large decrease of cosmic ray 
intensity was observed at Pioneer 11 and 10 during this period (Fillius and 
Axford, 1984). Several major magnetic storms and two large Forbush decreases 
were also recorded at the earth during the same period. For the second event, 
a significant disturbance in the outer heliosphere was observed by Pioneer 11 
at about 7 au and Pioneer 10 at about 16 au, manifested in the propagating 
shock waves and a large cosmic ray decrease, as well as two major geomagnetic 
disturbances, beginning on April 10 and 30, respectively. (Van Allen, 1979; 
McDonald et al. , 1981; Intriligator and Miller, 1982; Burlaga et al. , 1984). 

The purpose of this paper is to attempt to simulate qualitatively the 

* 

disturbances in the heliosphere during these two active periods using the 
* method developed by Hakamada and Akasofu (1982), providing a first order 
construction, temporally and spatially, of flare-generated shocks and their 
multiple interactions with each other, as well as with corotating interaction 
regions. Other dynamic, thermodynamic, and magnetic properties (other than 
first-order IMF distortion) are not simulated and can only be found from the 
MHD solutions. 

It is hoped that such a first order effort will be of some use in 
interpreting solar wind and cosmic ray observations by space probes. The 
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basis of this method is the distance (R) - time (t) relationship which can be 
obtained from the empirically constructed velocity (V) - distance (R) 
relationship from space probe observations. A theoretical V-R relationship 
(based on MHD solutions) can also be used to infer the R-t relationship. 
Olmsted and Akasofu (1984) demonstrated in detail the compatibility of our 
method with the MHD methods as far as the V-R and R-t relationships are 
concerned. Akasofu and Hakamada (1983) have already tested their method for 
successive six hypothetical flares from the same active region. 

2. Event of June - August, 1982 

Figure 1 shows the Hx synoptic charts for Carrington Rotations 1722, 1723 
and 1724 (Solar-Geophysical Data, no. 455, 456, 457 and 458). The three 
active regions are labeled as A, B and C. These regions were relatively short 
lived. Region A disintegrated after Carrington Rotation 1722, so that it is 
not seen in Carrington Rotation 1723. Region B was weakened during Carrington 
Rotation 1724. During Carrington Rotation 1724, Region C became the most 
active region. The Ha filtergram profiles around the central meridian passage 
of the three regions A, B and C are shown in Figures 2a (June 8), 2b (June 21) 
and 2c (June 18), respectively. Table 1 lists some of the major flares 
produced by the three regions and the storm sudden commencements (ssc's) which 

•i 

were inferred to be associated with the flares. 

Among the flares, we have chosen those which could, with reasonable 
confidence, be associated with an identifiable interplanetary shock wave at 
ISEE-3 and a ssc at the earth. The chosed flares are numbered, 1, 2, 3. ..7 as 
listed in Table 1; note that some of the flares are not numbered, because it 
was not possible to identify the resulting ssc for them; as we shall see 
later, it is not possible to simulate the associated event, since we cannot 
infer the flow speed from the transit time. However, the 'flares' no. 4a and 



4b are not observed ones and will be discussed later. Figure 3 shows the 
location of the flare, Venus (Pioneer 12), Earth, Pioneer 10 and 11 at the 
time of the eight flares. In this paper, our heliographic longitude is fixed 
in a heliospheric inertial frame. This frame is aligned with the commonly 
used heliocentric equatorial system, HEQ. In the HEQ system, the z axis is 
parallel to the solar rotational axis, x is directed along the intersection of 
the equatorial and ecliptic planes. This intersection is about 75° from the 
point Aries in the direction of the earth's orbital motion. More precisely 
HEQ x-axis is rotated by 74° 22' + 84' T, where T is the centuries since 
1900. The y-axis is formed by the right hand rule. Our heliographic 
longitude = 0° line (the x-axis) is rotated by an angle a with respect to the 
HEQ x-axis. Ue choose our heliographic longitude to coincide with Carrington 
longitude at the simulation start time, Tp « 0. Of course, as the simulation 
time progresses, the sun rotates and carries the Carrington coordinate system 
along, while our heliographic coordinates are fixed in space. As a point of 
reference, for the 1978 simulation, our heliographic longitude «* 0° line lies 
at HEQ longitude *» 101°, or heliocentric ecliptic, HEC, longitude => 176°, HEC 
latitude ** 7.1°. In the 1982 simulation our 0° line lies at HEQ longitude ■ 
162°, and HEC longitude 238®, HEC latitude => 2.2“ . 

A wide longitude and distance coverage of monitoring disturbances in the 
heliosphere make this particular period unique. Figure 4 shows, from the top, 
the solar wind speed observed at Pioneer 10, 11 and 12 (A. Barnes, private 
communication) and the magnetic field magnitude observed at ISEE-3 (E. J. 
Smith, private communication). One of the most remarkable features in Figure 
4 is the two shock waves observed at Pioneer 10 on July 30 and at Pioneer 11 
on August 3, respectively. Figure 5 shows the energetic particle data at 
Pioneer 10 and 11 and the neutron monitor record at Deep River for the period 
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between 1980 and 1983 (Fillius and Axford, 1984). Note in particular de- 
creases of the energetic particle fluxes at Pioneer 10 and Pioneer 11 in July 
and August, 1982 and the two large Forbush decreases at the earth in July. 
Recent new observations by deep space probes have also revealed that shock 
waves generated by successive solar activities advance well into the outer 
heliosphere and cause considerable disturbances of the solar wind and of the 
interplanetary magnetic field structure, resulting in reduction of the 
intensity of galactic cosmic rays (Van Allen, 1976; McDonald et al. , 1981; 
Smith, 1982; Burlaga, 1982; Burlaga et al. 1984). This figure will be 
discussed in detail after describing our simulation results. 

It is of great interest to follow how the shock waves generated by .the 
flare activity disturbed the heliosphere at two widely separated points, as 
well as at the earth. During the last several years, considerable progress 
has been made in understanding the propagation of solar wind disturbances in 
the heliosphere. The exploration of the outer heliosphere by deep space 
probes has recently added much new information on the general conditions of 
the solar wind and the magnetic field. In particular, the propagation of 
solar wind disturbances in the heliosphere has recently been studied 
extensively, both theoretically and observationally (cf. Dryer, 1974, 1975, 
1984; Dryer and Steinolfson, 1976; Smith et al. , 1977; Wu et al. , 1977, 1983; 
Dryer et al. , 1976, 1984; Intriligator, 1977; Smith and Wolfe, 1979; Wu, 1980; 

Wu et al. , 1977, 1979, 1983; Han et al. , 1984; Gislason et al. , 1984). 

The theoretical studies of the propagation of flare-generated shock waves 
are based on hydrodynamic or MHD methods which require considerable computa- 
tion. To simulate such complicated features in the outer heliosphere requires 
simultaneous solution of the MHD equations by finite-differencing techniques. 

At the present time, however, it is difficult to simulate by the MHD method an 
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extensive event associated with a large number of successive flares, superim- 
posed on the corotating interplanetary regions (CIR). It is for this reason 
that Hakamada and Akasofu (1982) devised a simple 'kinematic* method in 
simulating some aspects of the disturbed solar wind. 

During the course of our study, we have found that it is difficult to 
explain some of the major outer heliosphere events at Pioneer 12 on July '2 and 
at Pioneer 10 on July 30 only on the basis of the observed six flares, as well 
as of some less intense flares which were not listed in Table 1. We infer 
thus that there must have been intense flares behind the solar disk at Region 
A on about June 20 and at Region B on about July 1. There is no reason why 
active regions have to produce flares only when they are facing toward the 
earth. Further, there was no other major active region to cause heliospheric 
disturbances in the direction of Pioneer 12 and Pioneer 10 during this 
particular period. By including such hypothetical flares (no. 4a and 4b in 
Table 2), we attempt to explain the solar wind observations at Pioneers 10, 
and 12 and IMF observations at ISEE-3. 
a. Basic Flow Pattern 

We simulate first the basic flow pattern into which the disturbances 
generated by the seven flares are introduced. It is assumed that the 
heliomagnetic equator determined by Hoeksema et al. (1983) for the Carrington 
rotation 1722 (May 19-June 14) can be approximated by a two-stream (or so- 
called ' two-sector' ) situation and that the heliomagnetic equator remained the 
same during the Carrington rotations 1723 and 1724 (June 15-July 11 and July 
12-Augu8t 8); this is a reasonable assumption on the basis of inspection of 
the heliomagnetic equator for the rotations 1722, 1723 and 1724, determined by 
Hoeksema et al. (1983). In our simulation procedure, this situation is 
equivalent to assuming that the solar "dipole axis" is inclined with respect 
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to the rotation axis by angle x which is taken to be 20° in this particular 
period or that the heliomagnetic equator is given as x s ^ a (<t> + $ ) with 
respect to the heliographic equator where $ and $ o denote the longitude and 
phase angle, respectively. The solar wind speed V is assumed by Hakamada and 
Akasofu (1982) to increase towards higher magnetic latitude 6 as follows: 

V (km/sec) - 700 (1 - S 0 Vh\o\o 6 'f e 'f 7 + 300 

V (km/sec) =• 775 

The upper panel in Figure 6a shows the solar wind speed distribution on 
the source surface (of a spherical shell of 2.5 solar radii). In this way, 
the earth (located at latitudes 0 _<_ 7°) encounters a flow from the northern 
latitude during one-half of Carrington Rotation and a flow from the southern 
hemisphere during the other half. The flow pattern thus generated produces 
the familiar spiral interplanetary magnetic field (IMF) pattern with the 
alternating (toward and away) polarity, together with the two 'spiral arms', 
namely the CIRs. 

In this paper, we shall see how such a basic IMF pattern was disturbed by 
a series of solar flare-generated shock waves which took place in June and 
July 1982. We construct also the velocity (V) - time (t) pattern at Pioneer 
10, 11 and 12 and compare it with the observations. For details of the 
procedure of establishing the basic flow pattern, see Hakamada and Akasofu 
(1982). For the rest of the panel of Figure 6, see (b). 
b. Solar Wind Disturbances Induced by Flares 1, 2, 3, 5, 6 and 7 

In our simulation procedure, each flare is characterized by seven para- 
meters; they are the onset time Tp, heliographic longitude <p • latitude 0, the 
maximum speed Vp generated by the flare, the time variations of V character- 


for 0 < |0 | < 30° 
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ized by t, the extent characterized by the half width of the Gaussian distri- 
bution o° and the deceleration parameter a (Hakaoada and Akasofu, 1982). As 
an example, the middle panel in Figure 6a shows the solar wind speed distri- 
bution generated by the flare no. 1, which is characterized by $ « 302.4°, 6 » 
-12°, Vp - 820 km/sec, and a - 80° at the maximum epoch. The time variations 
of the maximum speed for the first flare (t - 5 hrs.) are shown at the bottom 
panel of Figure 6a; note that the speed reaches the maximum value C820 km/sec) 
five hours after flare onset. 

Table 2 gives the seven parameters chosen for our flares 1, 2, 3, 5, 6 
and 7, as well as for the hypothetical flares 4a and 4b (which will be dis- 
cussed in section c) . The time for our simulation study is reckoned from 1900 
UT, May 18, the zero-th day of Carrington Rotation 1722. For each flare, the 
onset time Tp, the longitude <j> and latitude 6 can directly be determined from 
the observations. The speed Vp, the deceleration parameter a, as well as T 
and o , must be determined on a trial and error basis in such a way that the 
arrival time of each shock wave at the earth and Pioneer 10 and 11 agrees 
approximately with the observations. Obviously, the choice of the parameters 
is not unique. In our method, there is no way, at the present time, to deter- 
mine Vp, t, a and a (or any other combinations of different sets of parameters 
which characterize a flare) uniquely. Any modeling effort of this kind would 
suffer from a similar non-unique problem with different degrees. One of the 
difficulties associated with this problem can be realized by noting that many 
flares occur far from the central meridian and thus that only, the skirt of the 
shock wave reaches the earth. The shock speed may depend greatly on the angle 
between the solar center-flare line and the direction of propagation, and this 
dependence may be different for different flares. Figure 6b shows the total 
(quiet plus flare-generated) solar wind speed distribution at the maximum 
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epoch of each of the eight flares Including that for flares no. 4a and 4b. 

Figures 7a-f show the solar wind disturbances caused in the inner 
heliosphere (< 2 au) by the flares 1, 2, 3, 5, 6 and 7 at the time of the 
arrival of the shock wave (at the time of the ssc for each flare) at the earth 
which Is indicated by the mark*. None of the six flares occurred near the 
central meridian of the solar disk, so that only the skirt of the shock passed 
by the earth. The shock waves must often extend well beyond ±45° in longitude 
from the flare longitude in order to explain the arrival at the earth (cf. 
Akasofu and Yoshida, 1967; Chao and Lepping, 1974; Dryer, 1974). 

Figures da and 8b show the disturbed solar wind in the outer heliosphere 
by the six flares; 1, 2, 3, 5, 6 and 7 at T ° 1752 hr (73.0 day, July 30) and 
T - 1848 hr (77.0 day, August 3), respectively (Figure 4). The location of 
Pioneer 10 ( ^ 28 au) and Pioneer 11 ( ~ 12 au) are indicated by the * mark. 
These dates correspond to the dates of the observation of the shock wave at 
Pioneer 10 and Pioneer 11, respectively; see Figure 4. One can see the pas- 
sage of a 'coalesced' shock at the location of Pioneer 11 at T ° 1848 hr (77.0 
day, August 3). However, it is not possible to explain the Pioneer 10 obser- 
vation at T =* 1752 hr (73.0 day, July 30) on the basis of the observed six 
flares, 1,2, 3, 5, 6 and 7. Region C caused a few other flares after July 14 
in the western hemisphere of the solar disk. In order for them to be respon- 
sible for the solar wind and energetic particle observations at Pioneer 10, 
unreasonably high speed flows are required, since the shock waves would have 
traveled about 28 au in less than 20 days, 
c. 'Flares 4a and 4b' 

One of the possiblilities to overcome this difficulty is to assume that 
there occurred two intense flares behind the visible disk of the sun on about 
June 20 and July 1, respectively. One of the reasons for suspecting it was 
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that the first active region A went behind the western limb on June 14, and 
the second active region B went behind the western limb of the sun on June 
27* From Figure 1, it is difficult to find any other active regions as the 
source(s) in explaining the intense shock waves observed at Pioneer 12 and 
Pioneer 10. Both Regions A and B were quite active near the western limb. 

The simulation study is complicated by the fact that no major increase of 
the "solar wind was observed on June 21 or 22 at Pioneer 12 which was located 
only 40° in longitude from Pioneer 10. It appears thus that the shock wave 
from Region A must have occurred on or before June 20 and must have been 
narrowly confined in longitude. We assume that a hypothetical flare designed 
as flare no. 4a occurred on June 16, although such as assumption is a very 
tentative one. 

There was also a large increase of the solar wind speed at Pioneer 12, 
from 469 km/sec on July 1 to 628 km/sec on July 2. Since the only active 
region facing Pioneer 12 was Region B, we infer that there occurred a fairly 
intense flare on July 1, which we designate as flare no. 4b. Another suppor- 
ting evidence of this hypothetical flare is the IMF disturbances observed by 
ISEE-3 on July 5-6 (see Figure 4). There were no significant solar activity 
between July 1 and July 5, so that it is possible that this particular distur- 
bance was caused by the hypothetical flare no. 4b; the skirt of the simulated 
shock wave must have passed by the earth on July 5. This situation is illus- 
trated in Figure 9a and 9b in which the location of the shock wave at T =» 1133 
hr (47.2 day, July 5) is shown. Unfortunately, it is not possible to associ- 
ate this particular flare as the cause of the event at Pioneer 10 on July 30, 
since an exceedinly high speed flow and the resulting shock wave are required. 

Figures 10 show the heliospheric disturbances in the ecliptic plane at T 
- 1437 hr (59.9 day, July 17), 1651 hr (68.8 day, July 26), 1752 hr (73.0 day. 
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July 30) and 1848 hr (77.0 day, August 3) by including the hypothetical 
flares 4a and 4b. Thus, what appears to be a singe major heliospheric event 
in Figures 4 and 5 may result from a number of flares,- including some behind 
the disk. We attempt to substantiate this claim a little more quantitatively 
in (d). 

d. ' The Velocity (V) - time (t) Relationship 

Figure 11 shows the V-t relationship at the location of Pioneer 10, 11 
and 12, respectively as well as at the location of the earth. They may be 
compared with the corresponding data in Figure 4. The V-t curves at Pioneer 
10, 11 and 12 have fair resemblance to the observed ones, as far as the two 
shock waves are concerned. However, much of the corotating structures seem to 
degenerate at the location of Pioneer 10 (28 au), since they are not very 
evident in Figure 4. On the other hand, some of them can be seen clearly at 
-the location of Pioneer 11 (13 au). It may be that many fine structures at 
Pioneer 12 are difficult to explain without including less intense solar 
activities. Note that it is also possible to identify tentatively some of the 
peaks observed at ISEE-3 with the computed peaks. 

3. Event of April - June. 1978 

Region 15266 of the sun caused a profound disturbance in the outer helio- 
sphere, and this subject has already been the source of several papers during 
the last several years. In particular, a flare of importance 3B at 1304 UT on 
April 28, 1978 was considered to be responsible for the intense shock wave and 
the large decrease of cosmic-ray intensity observed at Pioneer 11 at~7 au on 
May 11 and Pioneer 10 at~16 au on June 5 (Pyle et al., 1979; Van Allen, 1979; 
McDonald et al. , 1981; Intriligator and Miller, 1982; Burlaga et al. , 1984). 
Other intense shock waves were observed at Pioneer 10 on May 27 and at Pioneer 
11 on May 8. Actually, region 15266 had been quite active at least from the 
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day when It appeared near the eastern limb and caused a number of flares dur- 
ing the period when it was facing the earth (Dodson and Hedeman, 1981). There 
was at. least- another active region, region 15221, which caused several flares 
during the first week of April. 

Therefore, in order to study the major heliospheric disturbances in April 
and Jun,e 1978, it is important to examine how the shock wave generated by each 
flare propagated in the heliosphere and then how all the shock waves contribu- 
ted to the overall disturbances. Indeed, the intense Forbush decrease which 
began on April 30 showed a three-stage decrease; the main phase of the geomag- 
netic storm which began on April 30 had deepened until May 4. Both phenomena 
indicate that effects of the successive flares accumulated. 

In this paper we have chosen six intense flares of importance greater 

than 2 in April and May, 1978 and tentatively identified the corresponding 

ssc's; they are also listed in Table 3. Such identifiable ssc's are likely to 

be an important evidence for the generation of interplanetary shock waves. 

* 

Pyle et al. (1979) noted that the shock waves observed at Pioneer 10 on 
May 27 and at Pioneer 11 on May 8 cannot be explained without invoking an 
intense flare on. the backside of the solar disk on about April 15. Intriliga- 
tor.and Miller (1982) inferred that this particular flare occurred in region 
15266. The're was no highly active region within -±90° centered around region 
15266, so that this choice seems to be reasonable. However, this particular 
region' was at about 75° in heliographic (or Carrington) longitude, while the 
earth was at about 285°, so that in their Figure 1, l should be 27°, instead 
of 60°. This particular flare is designated as flare no. 3 in Table 3. How- 
ever, although it is very tentative, it is our finding that the space probe 
observations appear to be explained better by supposing that the responsible 
flare occured on April 17, instead of April 15. After the completion of this 
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study, we have learned that the Voyager- 1 and -2 spacecraft observed a very 
Important kilometric radio Type II event by the Planetary Radioastronomy Ex- 
periment on April 17, 1978; both spacecraft were behind the solar disk at that 
time (Y. Leblanc, private communication, 1984). Figure 12 shows the location 
of seven flares and of Pioneers .10 and 11, as well as of the earth, at the 
time of each flare listed in Table 3. 

The basic solar wind flow pattern in April - June, 1978, was different 
from that in June - August, 1982 in that the "four-sector" structure prevailed 
in the April - June, 1978 period (Hoeksema et al. , 1983) , whiie the "two- 
sector" structure did in the June - August, 1982 period. In our method, the 
"four-sector" structure is simulated by assuming that the heliospheric magne- 
tic equator is given by x sin (2$ + <t> 0 )» instead of x sin (<|> + <J> Q ). Four high 
speed streams and four "spiral arms" result from such a situation, instead of 
two high speed streams and two "spiral arms" in the "two-sector situation. In 
this paper, x is taken to be 20°. 

We introduce the seven flares into the "four-sector" situation thus de- 
termined. Table 4 gives the seven parameters for each flare. In this paper, 
the deceleration parameter a is taken to be 1000 for all the flares except for 
no. 4 flare (a = 2000). The flare paramenters must be determined on a trial 
and error basis, and in particular the April 28 flare (no. 4) had a very 
large a value (=190°). Such a large value of <J is needed, if this particular 
flare was responsible for the shock wave observed at both Pioneer 11 on May 11 
and Pioneer 10 on June 4. 

Figure 13 shows the geometry of the six shock waves generated by flares 
nos. 1, 2, 4, 5, 6, and 7, respectively, and their reslationship to the four- 
sector structure at the time when they reached the earth (causing ssc's). 
Since flare no. 3 occurred almost on opposite side of the earth with respect 
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to the sun, it is difficult to identify the associated ssc. Intriligator and 
Miller (1982) suggested that an ssc on April 17 was caused by the flare no. 
3. However, such an identification should be considered only very tentative; 
there was some indication that weak recurrent geomagnetic storms occurred 
successively about 27 days apart during this period, beginning on January 18, 
and then February 25, March 22 and April .17. 

One of the interesting features in Figure 13 is that the inner helio- 
sphere was almost completely surrounded by the “coalesced" shock waves, namely 
the four spiral arms and the shock wave generated by flare no. 4 (April 28) 
and could be an important cause for the major Forbush decrease. Such a unique 
situation was caused partially by a large value of o for flare no. 4 and could 
be an important qause for the major Forbush decrease. 

Figure 14 shows “snapshots" of the propagating heliospheric disturbances 
on May 8, 11, 27 and June 5. Those are chosen on the basis of Table 1 in In- 
triligator and Miller (1982). It should be recalled that the shock waves were 
observed at Pioneers 10 on May 27 and June 5 and at Pioneer 11 on May 8 and 
11. In this simulation study, it is assumed that there was no flare activity 
after May 1. As a result, the four-sector pattern began to develop in the 
inner heliosphere after May 1. Actually, there were two other intense flares 
(May ll.and May 31) which disturbed the inner heliosphere. Thus, the present 
simulation may serve in illustrating how a quiet condition might reappear 
after the end of an intense solar activity. 

Figure 15 compares the observed velocity profiles at Pioneers 10 and 11 
with the computed ones. The observed velocity profiles are taken from Intril- 
igator and Miller (1982). The major features of the velocity variations are 
fairly well reproduced at both Pioneers 10 and 11, indicating that our simula- 
tion may be a fair representation of the actual conditions to a distance of 
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about 7-10 au. However, the shock wave generated by . the April 28 flare ap- 
peared to degenerate considerably at the distance of Pioneer 10 (16 au). 
Thus, the shock structure which passed Pioneer 10 on June 5-13 could not be 
reproduced, although the decay part (after June 13) may have some resemblance. 
4. Discussions and Conclusions 

It is hoped that the simulation method adopted here is of some use in 
providing some idea about the geometry of the heliospheric disturbances in 
June - August, 1982 and April - June, 1978. It is possible that the helio- 
sphere was much more disturbed by other unseen major flares and weaker flares 
which are not included in this study. For the June - August, 1982 event, the 
simulated velocity- time profile at both Pioneers 10 and 11 has a fair resem- 
blance to the observed ones, suggesting that the major disturbances are simu- 
lated with fair accuracy. For the April - June, 1978 event, the fair agree- 
ment between the observed and simulated velocity variations at Pioneer 11 sug- 
gests that the simulated heliospheric disturbance patterns in the outer helio- 
sphere to a distance of about 10 au are a fair first approximation; however, 
the shock appeared to degenerate considerably by the time it reached Pioneer 
10 (16 au) , so that our simulation results are likely to be much less accurate 
at distances greater than 15 au. 

The two solar events produced a profound effect in the cosmic ray inten- 
sity observed at Pioneers 10 and 11, as well as at the earth (Figure 5). The 
cosmic ray variations during the April - June, 1978 event was described in 
detail by Van Allen (1979) and McDonald et al. (1981). Extending the impulse 
response function analysis by Bowe and Hutton (1982), Akasofu et al. (1984) 
have recently shown that a high (monthly) occurrence of major flares maintains 
some identity, in terms of its effect on cosmic-ray intensity, for as long as 
17 months. One could envisage that the disturbed structure which has some 
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resemblance to what we have simulated in this paper is propagated outward in 
the outer heliosphere for more than one year. It appears that the accumulated 
effects of such solar activity may have an important contribution to the 11- 
year modulation in the cosmic-ray intensity. In fact, one of the motives for 
s'tudying the June - August, 1982 episode is the effect which it had on gal- 
actic cosmic rays. During this short interval the cosmic ray intensity de- 
creased by an amount equivalent to half the amplitude of the eleven year solar 
cycle (Fillius et al. , 1983; Fillius and Axford, 1984). Thl3 feature is well 
illustrated in Figure 5. We see two declining epochs and two recovery 
epochs. The 1980 portion of this graph, up to the first cosmic ray minimum in 
early 1981, is part of a prolonged decline which began in 1978, from solar 
minimum conditions to solar maximum (Burlaga et al, 1984). The apparent 
recovery phase beginning in 1981 brought the intensities about halfway back to 
solar minimum values before being interrupted by the June - August, 1982 
events. The decrease in this short interval took the cosmic ray intensities 
back down to minimum values. Finally we see the more enduring increase that 
has lasted to the present time. Solar particle events interrupted the Pioneer 
11 curve in May, 1981, and again from June 6, 1982 until the decrease in early 
July.' 

■The first decline appears to be the sum of many stepwise events of dif- 
ferent magnitudes, which propogated outward at approximately the solar wind 
velocity, and frequently correlated with solar wind features (McDonald et al. , 
1981; McKibben, et al., 1982; Webber and Lockwood, 1981). However, the nature 
of these features remains elusive, with many suggestions such as bubbles, 
shocks, corotating interaction regions, and transient interaction regions 
(Newkirk et al. , 1981; Perko and Fisk, 1983; Burlaga et al. , 1984). 

The second decline is the largest concentrated change in at least a de- 
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cade siace the previous solar maximum in 1972. ' From the positions of the 
Pioneer spacecraft (Figure 3), we can see that the effect is large, not Just 
in amplitude, but in spatial extent. The Pioneer spacecraft are on opposite 
sides of the sun, and the decrease . evidently engulfed the entire helio- 
sphere. It is evident that we need to look at solar wind structures on the 
same scale, and in three dimension, to visualize the processes affecting the 
cosmic rays. The kinematic method used in this paper is, at present, the only 
method available to model solar wind interaction regions on this scale. 
Whether or not these images are accurate in all detail, they do demonstrate 
that, when many flares occur during one or more solar rotations, they can 
produce interplanetary structures on a global scale. It is quite plausible 
that such global structues have an especially strong cosmic ray modulation 
effect. 

Table 5 lists the times of occurence of some of the features that are 
identifiable in the data (Figures 4 and 5) and in the model (Figures 8, 10 and 
11). Of course, the model matches the solar wind data, by design, in as much 
detail as the modeling has been carried out. At certain times, there is a 
timing mismatch between solar wind and cosmic ray features that are certainly 
related (namely, the August events at Pioneer 11). This is not a modeling 
defect, because it is visible in the very data. It would be of Interest to 
carry out a detailed study of the solar wind structures in high time reso- 
lution and to try to clarify, using magnetic field data, the relationship 
between the solar wind and cosmic ray features. However, such a study is 
beyond the scope of this paper. 

Also, it is quite noticeable that, followng the initial Pioneer 10 cosmic 
ray decrease of July 30 - August 3, there are two more step decreases approx- 
imately 25 and 50 days later. These echo decreases do not seem to correspond 
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to any features in the solar wind veLocity at Pioneer 10 - another interesting 


puzzle. In the absence of local structures to produce the modulation, we 
should be forced to hypothesize non-local structures, for which there are 
certainly not enough data to constrain our models. Unfortunately, there are 
no data from the Pioneer 10 magnetometer to aid our understanding. The most 
likely explanation seems to be that the cosmic rays are modulated by magnetic 
field structures not evident in the solar wind velocity data. 
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Table 3 


McMa th 
Plage 



. • Date 

Flare 

Location 

M umber 

Importance 


88C 

1 

April 8, 1978 

0109 

M19 Wll 

15221 

28 

April 10 

1306 OT 

2 

April 11, 1978 

1334 

N22 W56 

15221 

2B 

April 13 

1925 

3 

April 17 

1200 

(unobserv) 

15266 




4 

April 28 

1304 

M2 2 E38 

15266 

3B 

April 30 

1925 

5 

April 29 

2010 

N20 E14 

15266 

2B 

May 1 

0828 

6 

4pril 30 

1420 

M2 8 El 4 

152166 

3B 

May 1 

1835 

7 

May 1 

1910 

N21 W12 

15266 

2B 

May 2 

2318 



Table 4 



Date 

flr(t) 

Longitude 

Latitude V 

T 

0 

t f 

1 

April 8 

0.0 

32.7 

19.0 

450.0 , 

10.0 

30.0 

1000.0 

2 

April 11 

84.4 

80.7 ' 

22.0 

800.0 

10.0 

60.0 

1000.0 

3 

April 15 

(226.8 

210.6 

22.0) 

400.0 

10.0 

95.0 

1000.0 

4 

April 28 

491.9 

3.3 

22.0 

580.0 

20.0 

190.0 

1000.0 

5 

April 29 

523.0 

28.3 

20.0 

900.0 

10.0 

60.0 

1000.0 

6 

April 30 

541.2 

29.2 

28.0 

1300.0 

10.0 

60.0 

1000.0 

7 

May 1 

570.0 

56.2 

21.0 

1100.0 

10.0 

90.0 

1000.0 



Table 5 


Time 

6/6-7 

6/24-28 

7/1 

7/2-3 

7/3-12 

7/6 

7/7-10 

7/9-10 

7/20-22 

7/25-26 

7/26 

7/29-30 

7/30-8/3 

8/1-3 

8/3 


Event at Pioneer 10 


CIR (model) 


CIR (model) 


CIR (model) 


Transient from flare 4a 
Step increase in Vsw 

Cosmic ray decrease 

CIR (model) 


8/7-9 


Event at Pioneer 11 
Solar Protons 


Vsw transient from 
flares 1 (model) 

Small increase in Vsw 

Cosmic ray decrease 

Vsw transient from 
flares 2 (model) 


Small increase in Vsw 


Small increase in Vsw 
CIR arrival (model) 


Vsw transient from 
flares 5 & 6 (model) 
Step increase in Vsw 

Cosmic ray decrease 


8/14-17 


CIR (model) 


Figure Captions 


Figure 1. 


Figure 2. 


Figure 3. 


Figure 4. 


Figure 5. 


Figure 6e. 


Hydrogen-Alpha synoptic charts £or Carrington Rotations 1722, 
1723 and 1724 from Solar-Geophysical Data, no. 455, 456, 457 
and 45S. The three active regions. A, B and C are indicated. 
Hydrogen-Alpha filtergram profiles near central meridian pas- 
sage time of ‘the three regions A, B and C (Solar-Geophysical 
Data, no. 455, 456, 457, 458). 

Satellite situation charts at the time of the 8 flares. The 
locations of the flare, Venus (Pioneer 12), Earth, Pioneers 

10, and 11 are shown at 00UT on each flare date. Note that in 
this figure only, the radial distance is given on a logar- 
ithmic scale in astronomical units (au) (A. Barnes, private 
communication). 

Solar speed observed at Pioneer 10, 11 and 12 and the magnetic 
field magnitude observed at ISEE-3 for the period 20 June to 
20 August, 1982 (E. J. Smith, private communication). 

Energetic particle Cherenkov Detector data at Pioneer 10 and 

11, and the Deep River neutron monitor record during the 
period between 1980 and 1983 (W. Fillius and V. 1. Axford, 
1984). 

Frame 1 and 2 show the background velocity distribution in 
heliographic coordinates (fixed on the source surface of 2.5 
solar radii) and the contribution to the velocity distribution 
due to flare #1, respectively. The third frame shows the time 
dependence of the particle velocity for flare #1. See text 
for details. 
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Figure 6b. 


Figure 7a-f. 


Figure 8a-b. 


Figure 9a-b. 


Figure 10a-d. 


Figure 11. 


Figure 12. 


Solar wind speed distribution at maximum epoch for the eight 
flares in Table 2. 

Equatorial plane projections of solar wind disturbances in the 
inner heliosphere (to 2 au) caused by the six flares 1, 2, 3, 
5, 6, 7 at the time of shock arrival at the earth. 

Equatorial plane projections of solar wind disturbances in the 
outer heliosphere (to 30 au) caused by the six flares 1, 2, 3, 
5, 6, 7 at T =* 1752 hr. (73.0 day, July 30), and T ■ 1848 hr. 
(77.0 day, August 3), respectively. The location of Pioneer 
10 (~28 au) in Figure 8a and Pioneer 11 (~12 au) in Figure 8b 
are indicated by the * mark. 

Equatorial plane solar wind disturbances at T 53 1133 hr (47.2 
day, July 5) with flare 4b included. Plots are shown out to 2 
au and 5 au respectively. 

Equatorial plane solar wind disturbances at-T ™ 1437 hr (59.9 
day, July 17), 1651 hr (68.8 day, July 26), 1752 hr (73.0 day, 
July 30) and 1848 hr (77.0 day, August 3) with flares 4a and 
4b Included. 

Simulated solar wind bulk speed at Pioneer 10, 11 and 12 and 
IMF magnitude at ISEE-3 for the period 20 June to 20 August, 
1982. Compare this figure tes in Table 2. . 

Satellite situation charts at the time of the 7 flares. The 
locations of the flare, Earth, Pioneers 10 and 11 are shown at 
00UT on each flare date. Note that the radical distance is 
given on a logarithmic scale in astronomical units (au). 
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Figure 13. 


Figure 14. 


Figure '15. 


Equatorial plane projections of solar wind disturbance in the 
inner heliosphere (to 2 au) caused by the six flares, 1, 2, 4, 
5, 6, 7 at the time of shock arrival at the earth. 

Equatorial plane projections of solar wind disturbance in the 
outer heliosphere on May 8, 11, 27 and June 5. The locations 
of Pioneers 10 and 11 are indicated. 

Observed and simulated solar wind bulk speed at Pioneers 10 
and 11. The observed speed is taken from Intriligator and 


Miller (1982) 
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